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Abstract 

An LRS Bianchi type-V cosmological models representing a viscous 
fluid distribution with a time dependent cosmological term A is inves- 
tigated. To get a determinate solution, the viscosity coefficient of bulk 
viscous fluid is assumed to be a power function of mass density. It turns 
out that the cosmological term A(t) is a decreasing function of time, which 
is consistent with recent observations of type la supernovae. Various phys- 
ical and kinematic features of these models have also been explored. 

PACS number: 98.80.Es, 98.80.-k 
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1 Introduction 

Cosmological models representing the early stages of the Universe have been 
studied by several authors. An LRS (Locally Rotationally Symmetric) Binachi 
type-V spatially homogeneous space-time creates more interest due to its richer 
structure both physically and geometrically than the standard perfect fluid FRW 
models. An LRS Bianchi type-V universe is a simple generalization of the 
Robertson- Walker metric with negative curvature. Most cosmological models 
assume that the matter in the universe can be described by 'dust' (a pressure- 
less distribution) or at best a perfect fluid. However, bulk viscosity is expected 
to play an important role at certain stages of expanding universe [T] — [5] . It has 
been shown that bulk viscosity leads to inflationary like solution [J and acts 
like a negative energy field in an expanding universe [S]. Furthermore, there 
are several processes which are expected to give rise to viscous effects. These 
are the decoupling of neutrinos during the radiation era and the decoupling of 
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radiation and matter during the recombination era. Bulk viscosity is associ- 
ated with the Grand Unification Theories (GUT) phase transition and string 
creation. Thus, we should consider the presence of a material distribution other 
than a perfect fluid to have realistic cosmological models (see Gr0n [B] for a 
review on cosmological models with bulk viscosity) . A number of authors have 
discussed cosmological solutions with bulk viscosity in various context [7J — [5]. 

Models with a relic cosmological constant A have received considerable at- 
tention recently among researchers for various reasons (see Refs.|10j — [Pf] and 
references therein). Some of the recent discussions on the cosmological constant 
"problem" and consequence on cosmology with a time- varying cosmological con- 
stant by Ratra and Peebles [15], Dolgov [16] — [18] and Sahni and Starobinsky 
[19] have pointed out that in the absence of any interaction with matter or radi- 
ation, the cosmological constant remains a "constant" . However, in the presence 
of interactions with matter or radiation, a solution of Einstein equations and 
the assumed equation of covariant conservation of stress-energy with a time- 
varying A can be found. For these solutions, conservation of energy requires 
decrease in the energy density of the vacuum component to be compensated by 
a corresponding increase in the energy density of matter or radiation. Earlier 
researchers on this topic, are contained in Zeldovich j2D], Weinberg [TT] and 
Carroll, Press and Turner [21]. Recent observations by Perlmutter et al. [22] 
and Riess et al. [23J strongly favour a significant and positive value of A. Their 
finding arise from the study of more than 50 type la supernovae with redshifts 
in the range 0.10 < z < 0.83 and these suggest Friedmann models with negative 
pressure matter such as a cosmological constant (A), domain walls or cosmic 
strings (Vilcnkin [24], Garnavich et al. [25]) Recently, Carmeli and Kuzmenko 
[2"6] have shown that the cosmological relativistic theory (Behar and Carmeli 
[27] ) predicts the value for cosmological constant A = 1.934 x 10~ 35 s~ 2 . This 
value of "A" is in excellent agreement with the measurements recently obtained 
by the High-Z Supernova Team and Supernova Cosmological Project (Garnavich 
et al. 25J, Perlmutter et al. [22], Riess et al. [23], Schmidt et al. [28 ). The 
main conclusion of these observations is that the expansion of the universe is 
accelerating. 

Several ansatz have been proposed in which the A term decays with time 
(see Refs. Gasperini [29]|30], Berman [31], Freese et al. [14], Ozer and Taha [14], 
Peebles and Ratra [32], Chen and Hu [33], Abdussattar and Viswakarma [34] . 
Gariel and Le Denmat [35], Pradhan et al. [3S]). Of the special interest is the 
ansatz A oc S~ 2 (where S is the scale factor of the Robertson- Walker metric) 
by Chen and Wu [33], which has been considered/modified by several authors 
( Abdel-Rahaman [37], Carvalho et al. [H], Waga [35], Silveira and Waga [35] . 
Vishwakarma [4"0]). 

Recently Bali and Yadav [41] obtained an LRS Bianchi type-V viscous fluid 
cosmological models in general relativity. Motivated by the situations discussed 
above, in this paper, we focus upon the exact solutions of Einstein's field equa- 
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tions in presence of a bulk viscous fluid in an expanding universe. We do this 
by extending the work of Bali and Yadav [41] by including a time dependent 
cosmological term A in the field equations. We have also assumed the coefficient 
of bulk viscosity to be a power function of mass density. This paper is organized 
as follows. The metric and the field equations are presented in section 2. In 
section 3 we deal with the solution of the field equations in presence of viscous 
fluid. The sections 3.1 and 3.2 contain the two different cases and also con- 
tain some physical aspects of these models respectively. Section 4 describe two 
models under suitable transformations. Finally in section 5 concluding remarks 
have been given. 

2 The Metric and Field Euations 

We consider LRS Bianchi type-V metric in the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 e 2x {dy 2 + dz 2 ), (1) 
where A and B are functions of t alone. 



The Einstein's field equations (in gravitational units c = 1, G = 1) read as 

1 



-M + Agj = -SttTV, 



(2) 



where R\ is the Ricci tensor; R = g 1 ^ Rij is the Ricci scalar; and Tf is the stress 
energy-tensor in the presence of bulk stress given by 

Ti = (p + p)viV j + pgl — {v\. + v J A + v j v e vi-i + ViV v? e )r] 



£- f 7 ? ) v Usl +v i v 3 ). 



(3) 



Here p, p, r\ and £ are the energy density, isotropic pressure, coefficients of 
shear viscosity and bulk viscous coefficient respectively and v % the flow vector 
satisfying the relations 

gijvV = -1. (4) 

The semicolon (; ) indicates covariant differentiation. We choose the coordinates 
to be comoving, so that v l = 8\. 

The Einstein's field equations ^ for the line element |T]) has been set up as 



2i?44 1 

~B~ + B 27 ~A 2 ~ 



p — 2tj 



4i 

A 



A. 



A44 B44 A4B4 1 

^ + ^- + ^B--A^ = - 8n 



-A, 



(5) 
(6) 
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2AiB 4 B\ 3 

^- + ^-^ = -8-P-A, (7) 

The suffix 4 after the symbols A, £? denotes ordinary differentiation with respect 
to t and 

= 4 



3 Solutions of the Field Eqations 

In this section, we have revisited the solutions obtained by Bali and Yadav [UJ . 
Equations © - © are four independent equations in seven unknowns A, B, p, 
p, £, r) and A. For complete determinacy of the system, we need three extra 
conditions. 

Eq. ©, after integration, reduce to 

A = B k , (9) 

where A; is an integrating constant. Equations ([5]) and © lead to 

-B44 A44 B\ A4B4 ( B4 Aj\ 

~B~- — -Bl-^W=- 167rr] {-B-^)- (10) 

Using Eq. © in JTJJ), we obtain 

S + (^)/ = -, 

where B 4 = /(B). Eq. HTJ leads to 

where L is an integrating constant. Eq. (|12p again leads to 

S= (fc + 2)^T2 (fc 1 -fc 2 e" 16,r, '*) F ^, (13) 

where 

*> - is? < 14 > 



167T?7 ' 

N being constant of integration. From Eqs. © and (fl~3"|) . we obtain 



(15) 
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Hence the metric ([!} reduces to the form 

ds 2 = -dt 2 + {k + 2)^T2 (j^ _ toe- 16 "" 4 )^ dx 2 

+ e 2x (k + 2)^T3 (to - toe~ 167 " ,t ) 7 ^ 5 (dy 2 + dz 2 ). (17) 
The pressure and density of the model (fT7|) are obtained as 

(87r)(167TT7) toe" 16 ^* 



87rp : 



3(/fc + 2) 2 (to - toe" 16 ^*) 2 



to(fc + 2) 2 (4ry + 3£) - {fc 2 (4ry + 3£) 



-4fc(r/ + 30 + 2(5r, + 60}fc 2 e- 16 H + ^-A, (18) 

J [(fc + 2)(to - k 2 e- 16nr >t)]k+2 

(2jfc+l), N2 2 e -32n V t 

87Tp= -JkTW {W7rfl) k2 (to - Aae-i*"*)* 

+ " ^+ A - ( 19 ) 

[(fc + 2) (to - fee" 167 ^*)] fe + 2 



(20) 



The expansion 6 in the model ([17]) is obtained as 

_ (167rr ? )toe- 16Tr, ' t 
~~ (to - fc 2 e- 167r? ' t )' 



For complete determinacy of the system we have to consider three extra condi- 
tions. Firstly we assume that the coefficient of shear viscosity is constant, i.e., 
V = Vo ( sa y)- For the specification of A(t), we secondly assume that the fluid 
obeys an equation of state of the form 

P = 1P, (21) 

where 7(0 < 7 < 1) is a constant. 

Thirdly bulk viscosity (£) is assumed to be a simple power function of the 
energy density [42] — [45], 

C(0 = 6)P n ; (22) 

where £0 and n are constants. For small density, n may even be equal to unity 
as used in Murphy's work [JS] for simplicity. If n = 1, Eq. (|2"2"|) may correspond 
to a radiative fluid [47]. Near the big bang, < n < \ is a more appropriate 
assumption [35] to obtain realistic models. 

For simplicity and realistic models of physical importance, we consider the 
following two cases (n = 0, 1): 
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3.1 Model I: 



Solution for n = 



When n = 0, Eq. reduces to £ = £ = constant. Hence, in this case Eqs. 
fig]) and 1(15)1, with the use of (EI]), lead to 



8tt(1 + 7)/9: 



8ttM r 



3iV 2 



h{k + 2) 2 (4?7o + 3Co) - {fc 2 (4r? + 3£ ) 
(2fc+l)Af 2 4 



iV 2 



+ 4fc(r ?0 + 3£o) + 2(5?? + eCo)}*^- 16 '"' * 

Eliminating p(t) between Eqs. p^|) and (|23[) . we obtain 
8ttM 



7V^ + 2 



(23) 



(1 + 7)A = 



3iV 2 



+ 2) 2 (4 % + 3£ ) - {k\A m + 3£ ) 



+ 4k( m + 3Co) + 2(5^ + 6eo)}fc 2 e- 16 ^°*J + (2fc + 1) 7 — 



(1 - 3 7 ) 

2k 



, (24) 



where 



3.2 Model II: 



A/ = 167rfc 2 ?7oe- 167r ' ,ot , 
JV= (fc + 2)(fci - fc 2 e- 167r?)ot ), 
P= 2fc 2 + 2fc + 5, 
Q = k 2 + 4fc + 4. 

Solution for n = 1 



(25) 



When n = 1, Eq. (f2"2")l reduces to £ = £oP ■ Hence, in this case Eqs. (|18p and 
ifTQll. with the use of ll2il. leads to 



8tt P - 



16wM{2k 1 (k + 2) 2 T] - Pk 2 r)Qt 



~167rr/o^ 



} 



3 [(1 + 7 )iV 2 - M{h{k + 2) 2 £o - Qfc 2 ^ e- 16 ^*}] 
4iV^T2 _ (2k + 1)M 2 



(26) 



[(1 + 7 )iV 2 - M + 2) 2 £o - QfeCoe- 16 -""'}] ' 
Eliminating p(t) between Eqs. (fHH) and (|26[) . we get 

A = 16,M[2 kl (k + 2) 2 Vo - Pk 2m e-^} + + ; (1 =^4 + 

(1 + 7) N 2 (i + 7 )jv*+3 



M[fc 1 (fc + 2) 2 e -Qfc 2 eoe 



-167r?7oi" 



{4AT-E+2 _ (2k + 1)M 2 } 



(1 + j)N 2 [(1 + 7 )7V 2 - M{k x {k + 2) 2 ^ - Qfeeoe- 16 ^'*}] (2?) 

From Eqs. (|23|) and (|26|) . we note that is a decreasing function of time 
and p > for all time in both models. The behaviour of the universe in these 
models will be determined by the cosmological term A; this term has the same 
effect as a uniform mass density p e // = — A/47rG, which is constant in space 



G 



and time. A positive value of A corresponds to a negative effective mass density 
(repulsion). Hence, we expect that in the universe with a positive value of A, 
the expansion will tend to accelerate; whereas in the universe with negative 
value of A, the expansion will slow down, stop and reverse. From Eqs. (j2"4")l and 
(|2"T)) . we observe that the cosmological term A in both models is a decreasing 
function of time and it approaches a small positive value as time increase more 
and more. This is a good agreement with recent observations of supernovae la 
(Garnavich et al. [25] , Perlmutter et al. 22J, Riess et al. [23] , Schmidt et al. 

EE]). 



The shear a in the model (jTTJ) is given by 

_ (k-l)M 

The non-vanishing components of conformal curvature tensor are given by 



(28) 



C 2323 = -Ci4i4 = (fc 3A ? 2 )M [ fcM " 167T77 fcl(fc + 2)], (29) 



C 13 i 3 = -0> 42 4 = ^ 3 ^] M [167rr/ fci (k + 2) - kM] , (30) 



daia = -C3434 = (fc 3jV 1 2 )M [16^ fc 1 (fc + 2) - kM]. (31) 



Equations $2U§ and (US]) lead to 

a (k-1) 



V3(/c + 2) 



= constant. (32) 



The model (fP7j) is expanding, non-rotating and shearing. Since ? = conatant, 
hence the model does not approach isotropy. The space-time (|17[) is Petrov type 
D in presence of viscosity. 

4 Other Models 

After using the transformation 

ki - k 2 e- W7nlt = sin (16tt7/t), k + 2 = 1/16tt??, (33) 
the metric (JTTJ) reduces to 



ds 2 = - 



cos (167t?7t) 



-1 2 


sin (167t?7t) 




167T?7 



A:i — sin (167t?7t) 

sin (167T77T)" 1 '''""'' J 



2(l-327r77) 



dx 



2 



+ e 2x 



16777] 



(dy 2 + dz 2 ). (34) 
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The pressure (p), density (p) and the expansion (0) of the model are ob- 
tained as 



8717? : 



(167rr/) 2 {fci — sin (1677777") 



3 sin 2 (I677777") 

(167777) (877£){fci — sin (1677777-)} 
sin (167t?7t) 



2fci-2(l-487r?7+11527rV){fci-sin(167r?7r)} 



I6777/ 



sin (167r?7r) 



2(l-32ir?)) 



2(24777? - l)(167777) 3 {fci - sin (I677777-)} 2 n 

"^P = ■ 2 \ 1~ 

sin (1677777-) 



167T77 



sin (I6777/7-) 



4.1 Model I: 



(167T7]){fci — sin (1677777-)} 
sin (1677777-) 

Solution for n = 



-A, (35) 

2(l-327rr;) 

+ A, 
(36) 
(37) 



When n = 0, Eq. (1221) reduces to £ = £0 = constant. Hence, in this case Eqs. 
flgg) and (JMJ), with the use of ([21]), lead to 



2(167T77o) 2 M : 
877(1 +7)P= 77—2 



\ [h - P l M l ] 



(167rr/ )(87rgo)Mi 
3sin 2 (167T77oT) L ' sin(167777 7-) 

2(24^0 - l)(167r77o) 3 M 2 



sin (16777707-) 



+ 4Ni 

Eliminating p(t) between Eqs. (|3"6")) and (|38|). we obtain 

2(167T7 7o ) 2 Af L _ p , (167r ?7 o)(87rCo)Af 1 



(38) 



(1+7)A = 



3 sin 2 (1 677770 r) 11 ' ' sin(167T7/or) 
2 7 (247r77o-l)(167r77 ) 3 A/ 2 



+ (1 - 37)JVi + 



sin 2 (167r?7or) 



(39) 



4.2 Model II: 



Solution for n 



When n = 1, Eq. ([22)1 reduces to £ = £0/? ■ Hence, in this case Eqs. (|35[) and 
(1361). with the use of JELL lead to 



87rp 



2(167r77o) 2 M 1 [(fc 1 - P1M1) + 3(24t777 - 1)(167T77 )M 1 ] 
3sin(167r77oT)[(l + 7) sin(167r?7oT) — 1677770^0^1] 

47V! sin(167777 r) 
+ [(1 + 7) sin(1677?7or) - 16^^] ' 
Eliminating p(t) between Eqs. (j3"6")) and (|40j) . we obtain 

2(1677?7o) 2 Afi(fci - Pi Mi) 



(40) 



A = 



3sin(167777or)[(l +7) sin(167777oT) — 1677770^0 M\ 
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[3(16^oeo)Mi + (1 - 3 7 ) sin(167rr?or)] 
1 [(1 + 7) sin(16irr, T) - lemjD&Mi] 

2(24^? - l)(167r7 7o ) 3 Jkr i 2 [7(l + 7) sin(167r?7 r) - (1 - 7 )(16^ g )M 1 ] 
(1 + 7) sin 2 (167r7 ?0 T)[(l + 7) sin(167r?7oT) - (167rry e )A^i] 



(41) 




167T?7 



-1 2(1-32^7;) 



sin (lGTr^oT") 



i 3 ! = 1 - 487T77o + 1152tt 2 ?^. 
The shear (cr) in the model ([34]) is obtained as 



(42) 



a = 



(1 — 487r?7o)(167rr7o)[fci — sin(167rr?o''") 
V3 sin(167r?7o''") 



(43) 



The models descibed in cases 4.1 and 4.2 preserve the same properties as in the 
cases of 3.1 and 3.2. 

5 Conclusions 

We have obtained a new class of LRS Bianchi type-V cosmological models of 
the universe in presence of a viscous fluid distribution with a time dependent 
cosmological term A. We have revisited the solutions obtained by Bali and Ya- 
dav [41] and obtained new solutions which also generalize their work. 

The cosmological constant is a parameter describing the energy density of 
the vacuum (empty space) , and a potentially important contribution to the dy- 
namical history of the universe. The physical interpretation of the cosmological 
constant as vacuum energy is supported by the existence of the "zero point" 
energy predicted by quantum mechanics. In quantum mechanics, particle and 
antiparticle pairs are consistently being created out of the vacuum. Even though 
these particles exist for only a short amount of time before annihilating each 
other they do give the vacuum a non-zero potential energy. In general relativity, 
all forms of energy should gravitate, including the energy of vacuum, hence the 
cosmological constant. A negative cosmological constant adds to the attractive 
gravity of matter, therefore universes with a negative cosmological constant are 
invariably doomed to re-collapse 49J. A positive cosmological constant resists 
the attractive gravity of matter due to its negative pressure. For most universes, 
the positive cosmological constant eventually dominates over the attraction of 
matter and drives the universe to expand exponentially |50j . 

The cosmological constants in all models given in Sections 3.1 and 3.2 are 
decreasing functions of time and they all approach a small and positive value 
at late times which are supported by the results from recent type la supernova 
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observations recently obtained by the High-z Supernova Team and Supernova 
Cosmological Project (Garnavich et al. [25], Perlmutter et al. [22], Riess et al. 
[55] . Schmidt et al. [2"B]), Thus, with our approach, we obtain a physically rele- 
vant decay law for the cosmological term unlike other investigators where adhoc 
laws were used to arrive at a mathematical expressions for the decaying vacuum 
energy. Our derived models provide a good agreement with the observational 
results. We have derived value for the cosmological constant A and attempted 
to formulate a physical interpretation for it. 
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